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I.  INTRODUCTION 


A  Monte  Carlo  computer  program  EDPEC1  has  been  constructed  at 
the  Ballistic  Research  Laboratory  to  calculate,  as  a  function  of  mass  depth 
along  the  beam  axis,  the  local  energy  deposition  by  an  electron  beam  and 
its  associated  secondary  electrons  and  positrons.  Since  the  distributional 
accuracy  of  calculations  using  such  Monte  Carlo  procedures  varies  as  l />/7v , 
where  N  is  the  number  of  sample  cascades,2  efficient  calculational  pro¬ 
cedures  are  generally  needed  to  keep  the  calculational  times  at  an  accept¬ 
able  level  when  high  accuracies  are  required  of  the  answers. 

In  this  connection,  the  integral  equations, 

*T»TJ  =  l-§r  <•> 

and 

_  f  ■glTl  dT 

i-«  • <(  =  m»,i) ,  (2) 

are  frequently  solved  to  find  respectively  the  average  energy  ( T0  -  TJ  lost 
by  an  electron  or  positron  which  has  traveled  a  known  mass  distance  u  and 
the  kinetic  energy  7,  at  which  a  catastrophic  electronic  interaction  occurs. 
The  quantity  70  in  the  preceding  equations  is  the  initial  kinetic  energy  of 
an  electron,  7U  is  the  average  residual  kinetic  energy  of  an  electron  which 
has  traveled  a  mass  distance  «(70,7U),  and  R/V(0,1)  is  a  random  number 
picked  with  equal  probability  at  any  point  in  the  range  from  0  to  1.  A 
catastrophic  interaction  is  identified  here  as  either  a  radiative  transition  in 
which  a  bremsstrahlung  photon  with  energy  greater  than  0.1-MeV  is  emit¬ 
ted  or  an  electron-electron  scattering  in  which  a  knock-on  electron  having 
energy  greater  than  0.1-MeV  is  ejected. 


IW.B.  Beverly,  “EDPEC:  A  Computer  Program  for  Calculating  the  Energy  Deposition 
by  an  Photon-Electron  Cascade  in  Multi-Slab  Arrays,"  In  Publication  as  a  BRL 
Technical  Report. 

*Y.U.  Shreider,  The  Monte  Carlo  Method.  The  Method  of  Statistical  Trials.  Pergamon 
Press,  Aylesbury,  Bucks,  Great  Britain,  1967. 


The  quantity  e{T)  is  the  restricted  stopping  power*  (MeV-cm*/g)  of 
electrons  and  positrons  in  a  material,  and  o{T)  is  the  mass  cross  section 
(cm*/g)  for  bremsstrahlung  production.  All  functions  used  in  the  preceding 
equations  are  derived  by  piecewise  fitting  tabular  data  to  interpolation 
functions  tractable  for  calculations  of  the  type  described  in  this  report. 
Since  a  significant  fraction  of  the  total  calculational  effort  of  EDPEC  is 
expended  on  evaluating  these  integrals,  an  increase  in  the  associated  calcu¬ 
lational  efficiency  would  be  expected  to  significantly  reduce  the  overall  cal¬ 
culational  effort. 

These  integral  equations  were  initially  solved  using  the  bisection 
method4  where  about  twenty  iterations  were  generally  needed  to  reduce  the 
interval  enclosing  an  answer  to  0.0001  percent  of  its  initial  width.  How¬ 
ever,  since  the  calculational  times  for  problems  of  interest  turned  out  to  be 
very  long,  an  attempt  was  made  to  find  more  efficient  procedures.  VVe 
describe  such  procedures,  which  are  derived  by  applying  Newton’s  Method4 
for  finding  the  zero  of  an  analytic  function,  in  the  next  section  of  this 
report.  Example  calculations,  where  tractable  analytic  functions  are  substi¬ 
tuted  for  e( T)  in  Equation  1  and  <r{T)/e(T)  in  Equation  2,  are  described  in 
Section  III. 


#M.J.  Berger  and  S.H.  Seltzer,  ‘‘Stopping  Powers  and  Ranges  of  Electrons  and  Positrons 
(2nd  Ed.),”  Prepared  for:  U.S.  Department  of  Commerce,  Office  of  Health  and 
Environmental  Research,  Office  of  Naval  Research.  NBSIR  82-2550-A,  December  1982. 

4S.D.  Conte  and  C.  deBoor,  ELEMENTARY  NUMERICAL  ANALYSIS.  An  Algorithmic 
Approach.  McGraw-Hill  Book  Company,  Third  Edition  (1980). 


II.  THE  SOLUTION  PROCEDURES 


An  improved  estimate  xIi+1  of  the  location  xt  of  a  zero  of  the  equation 


F\x)  =  0 

can  be  derived  from  an  estimate  xj  by  using  Newton’s  Method4,  that  is 


- 


F'ti) 


(3) 

(4) 


if  the  first  derivative  F'(xj)  is  continuous  and  non-zero  between  x*  and  xt. 
The  application  of  this  method  to  the  aforementioned  EDPEC  integral 
equations  will  now  be  described. 


The  electron-range  problem  as  developed  in  EDPEC  is  defined  in  the 
following  manner: 


1.  An  electron  having  an  initial  kinetic  energy  T0  travels  a  known  mass 
distance  u(T0,Ta)  while  constantly  losing  energy.  The  unknown  residual 
kinetic  energy,  known  to  lie  between  T0  and  TXl  is  obtained  by  solving 
Equation  1, 

=  1&-  (,) 


The  range  of  validity  of  a  fitted  interpolation  function  extends  from  T0 
down  to  r,.  Since  the  form  of  t(I)  used  by  EDPEC  does  not  permit  the 


dT 


indefinite  integral  t°  be  expressed  in  a  closed  form,  any  associated 

definite  integral  is  evaluated  by  expanding  the  integrand  in  a  convergent 


dT 

series.  Consequently,  the  definite  integral  /  -  and  any  intermediate 

r, 

definite  integral  whose  lower  limit  lies  between  T0  and  r,  can  be  approxi¬ 
mated  to  any  specified  accuracy  by  evaluating  a  sufficient  number  of  terms. 


2.  The  electron  mass  range  integral  (Equation  1)  can  be  rearranged  to 
the  form  of  Equation  3,  that  is 

t, 

nTJ  =  l  =  0  (5A> 

The  first  derivative,  needed  in  the  Newtonian  iteration  formula,  is 


since  T0  and  v(T0,TJ  are  known  in  the  problems  being  discussed  and  Tu  is 
the  unknown  quantity. 

3.  The  derivative  of  T)  is  non-zero  on  the  range  from  T0  to  Tx  since 
i{T)  is  positive  and  non-zero.  Consequently,  Newton’s  Method  may  be 
applied  to  calculate  a  better  approximation  Tj*1  from  any  estimate  Tj 
which  lies  between  T0  and  Tx. 

T'  AT 

4.  The  mass  distance,  «(^0)^i)==/  35  we^  ^  the  associated 

T, 

indefinite  integrals,  have  been  calculated  earlier,  and  these  values  can  be 
retrieved  and  used  in  further  calculations  such  as  the  derivation  of  an  ini¬ 
tial  estimate  Tj  of  7*u. 

5.  The  value  of  any  element  in  the  set  of  quantities, 
[77,  e(77),  fir,'),  F\Tj),  «(r0,7V)l ,  needed  during  the  calculation  of  a  first  or 
improved  estimate,  has  been  or  can  be  easily  calculated. 


Item  5  completes  the  definition  of  the  electron-range  integral-equation 
problem.  A  first  estimate  T*  is  derived  by  assuming  that  the  stopping 
power  between  T0  and  Tx  is  approximated  by  the  linear  relation  cf(  T)  whose 
slope  m  is  derived  from  the  points,  [r0,  c(r0)]  and  [7*,,  c(r,)J  (Figure  1),  that 
is 


m 


c(  r0)  -  €(  r,) 


(6A) 


The  function  (*(7)  is  then  derived  by  translating  the  points,  [r0c(r0)]  and 
[7’|,c( 7*1 )],  a  distance  Ac  so  that  the  mass  distance  wt( T*0,  T,),  calculated  using 
cf(7) ,  is  equal  to  u(T0,  Tx),  that  is 


«f(  T0,TX)  =  / 


=  / 


r,  *'(  t)  f,  Mr- r,)  +  <(r,)  +  At 


=  «<r„7i).  (6B) 


The  magnitude  of  the  translation  is  derived  to  be 

<( r,) «" “(T,T,)  -  m  ( 7-0  -  r,)  - .( r,) 

Ac  =  - -  r-sr? - 

l  -  em  u*  J»’T,1 

and  the  linear  equation  i\T)  is  given  by 


«*m  =  m(T-  r,)  +  .(r,)  + At  . 


(6C) 

(01)) 


Figure  1.  The  Linear  Approximation  of  f(7)  Used  To  Derive  a  First 
Estimate  of  T%  in  the  Mass-Range  Problem. 
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The  value  of  the  first  estimate  Tu*  of  the  residual  energy  of  an  electron 
which  has  traveled  a  known  mass  distance  v(T0,Tu)  is  then  calculated  by 
solving  the  integral  equation 


r. 


/ 

t: 


_ dT _ 

m(T-  r1)  +  «(r1)  + A£ 


«( T0,  TJ  . 


(6E) 


This  is  found  to  be 


r„‘  =  r,  -  -i-  [<(T,)  +  A(  - «—  “(T,'Ti1  ( m  ( r,  -  rf)  +  d  r,)  +  )  ] .  (6F) 

Finally,  the  first  estimate  just  calculated  is  iteratively  improved  using 
Equation  4  until  a  specified  accuracy  has  been  attained. 


The  first  estimate  as  derived  above  has  proved  in  practice  to  be 
sufficiently  accurate  so  that  subsequent  iterations  rapidly  converged  toward 
the  true  answer.  In  general,  three  or  four  iterations  are  sufficient  to  attain 
an  accuracy  well  within  the  overall  EDPEC  calculational  accuracy. 


The  catastrophic  sampling  problem  as  developed  in  EDPEC  is  defined 
in  the  following  manner: 


1.  An  electron  having  an  initial  kinetic  energy  T0  loses  energy  continu¬ 
ously  as  it  travels  through  a  medium  until  it  suffers  a  catastrophic  interac¬ 
tion  during  which  a  finite  amount  of  energy  is  instantaneously  lost.  The 
kinetic  e^  ~gy  of  an  electron  as  it  enters  into  a  catastrophic  interaction  is 
obtained  by  solving  Equation  2, 

_  J  fill  dx 

1  -  e  ’•  <|T)  =  0,1) ,  (2) 

where  Tt  is  known  to  lie  between  T0  and  Tv  The  mass  cross  sections  for 
catastrophic  interactions  are  also  piecewise  fitted  on  the  interval  from  T0  to 
T,  by  functions  <r(T)  which  are  tractable  to  sampling  calculations  of  the 
type  discussed  in  this  report.  As  described  earlier  in  Item  1  of  the  mass- 

range-problem  description,  the  integrals  associated  with  Equation  2  cannot 

Tt 

be  given  in  closed  form.  Consequently,  the  definite  integral  /  dT  and 

any  intermediate  definite  integral  whose  lower  limit  lies  between  T0  and  Tx 
are  evaluated  to  an  acceptable  accuracy  by  summing  a  sufficient  number  of 
terms  of  a  convergent  series  expansion. 
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2.  The  catastrophic  sampling  integral,  (Equation  2)  can  be  rearranged 
to  the  form  of  Equation  3,  that  is 

f\T.)  = 

-  f  ^IldT 

1  -  «  1 1  1  -  RN[0,1)  =  1  -  e_c(T*>T«)  _  RN(ofi)  =  0  .  (7A) 

The  first  derivative  F  { T #)  of  F\Tt),  needed  by  Newton’s  iteration  formula,  is 


since  T0  and  RN(0,l)  are  known  constants. 

3.  Since  the  functions,  <t{T)  and  ((T),  are  positive  and  non-zero  between 
Tq  and  Tlt  the  function  FI(T)  does  not  vanish.  Consequently,  Newton’s 
Method  may  be  applied  to  calculate  a  better  approximation  r,,+1  when  T,' 
lies  between  T0  and  Tx. 

4.  The  definite  integral,  f  dT  =  c(r0,7\) ,  along  with  the  associ- 

ated  indefinite  integrals  at  T0  and  Tt  ,  have  been  calculated  earlier  and  can 
be  retrieved  for  subsequent  calculations  such  as  the  derivation  of  an  initial 
estimate  T*  of  T%.  For  calculational  convenience,  g{T)  —  <r (T)/c(T)  in  the  fol¬ 
lowing  discussion. 

5.  The  value  of  anv  element  in  the  set,  [77,  g(T,%  F\T,%  F’(Tj),  f(r0,rB)] 
between  T0  and  Tx  has  been  or  can  be  easily  calculated  during  the  deriva¬ 
tion  of  an  initial  estimate  or  the  conduction  of  subsequent  iterations. 

Item  5  completes  the  definition  of  the  catastrophic-sampling  integral- 
equation  problem.  A  first  estimate  Tj  of  the  solution  of  Equation  7A  is 
derived  by  approximating  the  function  g(  T)  =  <r(  T)/t(  T)  by  the  more  tract¬ 
able  linear  relation  g\T).  One  point  on  g*(T)  is  constrained  to  pass  through 
the  point  (r0,y(r0)]  (Figure  2),  and  the  approximation  integral  is  set 

equal  to  the  true  integral,  that  is 

r*  T%  tr\ 

c\T0Jx)  =  /  g\T\dT  =  ciT0,Tx)  =  /  ^  dT .  (8A) 

r,  r,  1 1 

Performing  the  rotation  about  T0  as  described,  the  value  of  g  at  7,  is 


3Z3? 


derived  to  be 


2Mr0,ri)l  _ 

r(Ti)  =  — = - r - rtro)- 

1  o  ~ 

The  linear  approximation  of  y(7)  from  T0  to  Tx  becomes 

M  =  m{T-  T0)  +  <AT0) 

whose  slope  m  is 

ni^-gHW 


m  — 


Ta-  r, 


(8B) 


(8C) 


(8D) 


The  first  estimate  T*  of  Tt  is  now  derived  by  solving  the  integral  equa¬ 


tion 


t  C  (T«-T.1)  _ 


which  can  be  changed  to 

r. 

f  [m(T-  T0)  +  «(r0)l  dT  =  -  /n[l  -  i?iV(0,l)]  . 
V 

This  reduces  to  the  quadratic  equation 


(8E) 


(8F) 


(T.1)5  -  2  I  To  -  |  ( r,1)  +  T}  -  -1  /n[l  -  /f/V(0,l)|  -  *l(^)  r°  =  0  (8G) 

mm  m 


of  the  variable  T,  whose  solutions  are  the  well  known 

(72)  = 

where: 


2  a 


a 

6 


and 


-2 


r„- 


J*7o) 


m 


c 


r0*  -  -1  /nil  -  /?A/10,1))  - 
m 


2  x  g(r0)  X  T0 

m 


It  is  noted  without  explanation  that  Equation  8G  is  never  solved  in  EDPEC 
for  RN( 0,1)  =  1  so  that  the  infinity  of  ln( o)  never  occurs  in  practice.  The 
appropriate  solution  to  the  problem  at  hand  is  selected  from  the  two 


available  solutions  as  that  solution  which  lies  between  T0  and  Tv  This  esti¬ 
mate  is  iteratively  improved  using  Newton’s  Method  until  an  acceptable 
accuracy  has  been  attained. 

The  first  estimate  as  derived  above  has  proved  to  be  quite  close  to  the 
true  answer  for  a  test  problem  whose  solution  is  illustrated  later  in  Section 
III  of  this  report.  Four  or  five  subsequent  iterations  generally  proved 
sufficient  to  provide  an  acceptable  accuracy  when  this  first  estimate  was 
inserted  into  Newton’s  Procedure  (Equation  4). 
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III.  EXAMPLE  SOLUTIONS  OF  THE  TWO  INTEGRAL  EQUATIONS 


Tractable  functions  were  selected  in  turn  forT/e(7)  (Equation  1)  and 
a(T)/i(T)  —  g(T)  (Equation  2)  which  are  illustrative  of  the  behavior  of  these 
functions  over  the  electron  energy  ranges  of  interest.  Example  problems 
were  devised  using  these  equations  and  a  short  computer  program  was  writ¬ 
ten  to  solve  each  problem.  The  range  of  values  of  the  dependent  variable  T 
was  selected  for  calculational  convenience  and  does  not  correspond  to  the 
kinetic  energy  (MeV)  of  the  electrons  in  problems  of  interest. 


The  stopping  power  in  the  mass-range  problem  is  taken  to  be 

t(T)  =  l+0.ir2  .  (9A) 

The  initial  energy  of  an  electron  is  taken  to  be  T—  3  (Arbitrary  Units)  and 
solutions  are  sought  for  a  range  of  mass  distances  up  to  that  associated 
with  a  total  loss  of  kinetic  energy  (T=  0).  The  mass  range  integral  is 


and  the  associated  equation  to  be  solved  is 

™  -  °' 


(9B) 


(9C) 


where  T0  and  u (T0,TU)  are  known  and  ru  is  unknown.  The  derivative  F'(rj 
of  /\ru)  is 


F'(TU)  =  -±-\F\TJ)  = 


i  +  o.i  r* 


(9E) 


The  maximum  mass  distance  traveled  by  an  electron  in  the  specified  energy 
range  is  identified  as 


«(3,0)  =  / - , 

o  1  +  o.ir* 


(9F) 


and  can  be  easily  evaluated  in  closed  form. 


Solutions  of  Tm  associated  with  mass  depths  ranging  from  the  very 
small  to  those  close  to  the  maximum  value  of  «(3,0)  are  needed  in  order  to 
demonstrate  whether  the  developed  procedures  are  effective  over  the  total 
energy  range  of  a  problem.  In  this  connection,  solutions  are  calculated  in 


turn  for  depths  given  by 


«n[3,(ru)n) 


(n-1)  43,0) 

25 


where  n  ranges  in  value  from  2  to  25  as  well  as  the  values, 


and 


tqfrfrj,]  =  0.0001  , 


ru),#]  =  «(3,0)  -  o.oooi . 


A  short  computer  program  was  constructed  to  solve  this  set  of  prob¬ 
lems.  The  first  step  in  calculating  a  solution  for  a  specified  mass  depth  is 
to  derive  the  linear  approximation  <*(T)  of  t(T)  (Figure  3).  The  slope  m  and 
translation  increment  A«  are  derived  as  described  by  Equations  6A  and  6D, 
respectively.  The  initial  estimate  T,1  is  then  calculated  using  Equation  6F. 
Finally,  the  initial  estimate  is  iteratively  improved  using  Newton’s  Method 
(Equation  4)  until  a  specified  accuracy  has  been  attained. 

Representative  results  from  the  example  calculation  are  given  in  Table 
1.  Convergence  is  very  rapid  for  each  problem  in  the  set;  three  iterations 
after  the  first  estimate  are  sufficient  in  all  cases  to  obtain  a  very  good  accu¬ 
racy. 


The  integrand  g(T)  associated  with  the  catastrophic  sainpling  problem 
is  taken  to  be 


1(71  =  7m  =  T’-3T+ 4  <10A) 

The  initial  energy  of  an  electron  is  taken  to  be  3  (Arbitrary  Units),  and  the 
transport  of  the  electron  is  to  be  conducted  until  it  either  suffers  a  catas¬ 
trophic  interaction  or  its  energy  is  reduced  to  1.  The  catastrophic  sampling 
equation  is 

l_e-c(T„T.)  =  RN[Ql)  (JOB) 

where 


c(T0,T.) 


r,  r, 

/  p(D  AT  =  /  (7*  -  3T+  4)  iT. 

t.  t. 


The  associated  equation  to  be  solved  is 

f\  r.)  =  i  -  rn[  o,i)  -  =  o 


(10C) 


( 10D) 


Figure  3.  The  Linear  Approximation  of  f(  T)  Applied  to  the  Example 
Mass-Range  Problem. 


u(3,0)  =  2.4003907662633  A«  = 

-0.1466557965344 

D 

T.1 

fitJi 

ATu1 

n  *  1 

u(3,T  )  =  0.0000100000000 

i 

2.9999824665843 

-0.0000007718609 

0.0000014665276 

2 

2.99998100U0567 

0.0000000000002 

-0.0000000000003 

3 

2.9999810000570 

-0.0000000000000 

0.0000000000000 

n  =  5 

u(3,T  )  =  0.3840625224421 

1 

2.363952U227107 

-0.0136311085117 

0.0212485388507 

2 

2.3427034838600 

0.0000439806607 

-0.0000681183890 

3 

2.3427716022490 

0.0000000004531 

-0.0000000007019 

4 

2.3427716029509 

-0.0000000000000 

0.0000000000000 

n  =  10 

u(3,T  )  =  0.8641406754948 

1 

1.6653316461874 

0.0034329650556 

-0.0043850393790 

2 

1.6697166855663 

0.0000019628597 

-0.0000025100959 

3 

1.6697191956623 

0.0000000000006 

-0.0000000000008 

4 

1.6697191956631 

-0.0000000000000 

0.0000000000000 

n  *  15 

u(3,T  )  =  1.3442188285475 

1 

1.0604190045084 

0.0330169617933 

-0.0367296810623 

2 

1.0971486855707 

0.0001163822100 

-0.0001303915468 

3 

1.0972790771175 

0.0000000014861 

-0.0000000016650 

4 

1.0972790787825 

0.0000000000000 

0.0000000000000 

n  *  20 

u(3,T  )  *  1.8242969816001 

1 

0.5366448399415 

0.0445132895760 

-0.0457952174943 

2 

0.5824400574357 

0.0001089978658 

-0.0001126954704 

3 

0.5825527529061 

0.0000000006920 

-0.0000000007155 

4 

0.5825527536216 

0.0000000000000 

0.0000000000000 

n  *  26 

u(3,T  )  *  2.4003807652633  1 

1 

0.0000085334648 

0.0000014665452 

-0.0000014665452 

2 

0.0000100000000 

0.0000000000000 

-0.0000000000000 

3 

0.0000100000000 

0.0000000000000 

0.0000000000000 

Table  1.  Representative  Solutions  of  the  Mass-Range  Example  Problem. 


and  its  derivative  F'(T%)  is 

fir.)  =  jrf\T.)  =  [-4^,7.)]  .  (10E) 


Since  the  maximum  value  of  the  catastrophic  sampling  integral  on  this 
energy  range  is 

3 

c(3,l)  =  /  (7«-3r+  4)  AT  =  14/3,  (10F) 

l 

then  a  catastrophic  interaction  will  only  occur  when  the  random  number 
lies  on  the  range 

0  <  RN{0,1)  <  1  -  «  l4/3  .  (10G) 

Solutions  of  T,  associated  with  random  numbers  ranging  from  those  close  to 
0  to  those  values  close  to  the  maximum  value  described  above  in  Equation 
10G  are  used  in  order  to  demonstrate  if  the  procedures  are  effective  over 
the  total  energy  range  of  a  problem.  In  this  connection,  the  random 
numbers,  ordinarily  selected  using  stochastic  procedures,  are  calculated 
deterministically  as 

RNn( 0,1)  =  ~  *)■  ~  ?~14-3)- 

25 

where  n  ranges  in  value  from  2  to  25  as  well  as  the  values 


and  for 


RNt(0,l)  =  0.0001 , 


RN2s{ 0,1)  =  1  -  «'14/3  -  0.0001  . 


A  short  computer  program  was  constructed  to  solve  this  set  of  prob¬ 
lems.  The  first  step  in  calculating  a  solution  for  a  specified  random  number 
is  to  derive  the  linear  approximation  f(1)  of  y(T)  (Figure  4).  The  derivation 
of  the  linear  approximation  used  here  differs  from  that  described  earlier  for 
the  mass-range  problem  in  that  q\T\  is  derived  by  rotating  u(T)  about  T0 
until  the  associated  sampling  integral  t»t(r0,7’1 )  is  equal  to  «(r0,r,)  (Figure  4). 
This  rotation  is  conducted  as  described  by  Equation  8B,  and  an  initial  esti¬ 
mate  Ttl  is  calculated  by  solving  Equation  8G.  Finally,  the  initial  estimate 
is  iteratively  improved  using  Newton’s  Method  until  a  specified  accuracy  is 
obtained. 
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Figure  4.  The  Linear  Approximation  of  <t(7V<(7’)  =  AT)  Applied  to  the 
Example  Catastrophic  Sampling  Problem. 
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Hiiiwrrirv  "v 


Representative  results  are  given  in  Table  2.  Convergence  is  very  rapid 
for  all  cases;  five  iterations  from  initial  estimate  are  sufficient  in  all  cases  to 
obtain  a  very  good  accuracy.  -  — 


c(3,l)  *  4.6666666666667  [g+(T)]' 

=  0.6666666666667 

I 

Ts: 

F(TSJ) 

ATs1 

n  =  5 

rn  =  0.1584952699918 

1 

2.9564642213854 

-0.0010408024021 

0.0003190945754 

2 

2.9561451268099 

-0.0000007673761 

0.0000002356134 

3 

2.9561448911965 

-0.0000000000004 

0.0000000000001 

4 

2.9561448911964 

-0.0000000000000 

0.0000000000000 

n  =  20 

rn  =  0.7528525324609 

1 

2.62U5632622800 

-0.0199883563891 

0.0248945802238 

2 

2.5956686820562 

-0.0009004277777 

0.0012303231939 

3 

2.5944383588623 

-0.0000020420717 

0.0000028029407 

4 

2.5944355559216 

-0.0000000000106 

0.0000000000145 

5 

2.5944355559072 

0.0000000000000 

-0.0000000000000 

n  =  24 

rn  =  0.9113478024526 

1 

2.2888927291901 

-0.0215122039163 

0.0823122638492 

2 

2.2065804653409 

-0.0024388256425 

0.0119032702411 

3 

2.1946771950998 

-0.0000411876600 

0.0002080029803 

4 

2.1944691921195 

-0.0000000122261 

0.0000000617803 

5 

2.1944691303393 

-0.0000000000000 

0.0000000000000 

6 

2.1944691303393 

0.0000000000000 

-0.0000000000000 

n  =  25 

rn  *  0.9509716199506 

1 

2.0634113718903 

-0.0178810679518 

0.1292631377944 

2 

1.9341482340959 

-0.0026371822739 

0.0263315613811 

3 

1.9078166727148 

-0.0000806562499 

0.0008570572678 

4 

1.9069596154470 

-0.0000000808750 

0.0000008611080 

5 

1.9069587543390 

-0.0000000000001 

0.0000000000009 

6 

1.9069587543381 

-0.0000000000000 

0.0000000000000 

n  »  26 

rn  =  0.9905854374485 

1 

1.0017498014171 

-0.0000219519556 

0.0011641557120 

2 

1.0005856457051 

-0.0000000191512 

0.0000010173989 

3 

1.0005846283061 

-0.0000000000000 

0.0000000000008 

4 

1.0005846283054 

-0.0000000000000 

0.0000000000000 

1 

Table  2.  Representative  Solutions  of  the  Catastrophic  Sampling  Problem. 
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IV.  CONCLUSIONS  AND  PROGNOSIS 


Improved  procedures  for  solving  certain  integral  equations  associated 
with  the  Monte  Carlo  simulation  of  electron-photon  cascades  have  been 
investigated.  The  efficiency  of  these  procedures  has  been  illustrated  for 
integral  equations  where  tractable  functions,  representative  of  the 
bremsstrahlung  cross  sections  and  the  electron  stopping  power,  have  been 
substituted.  In  fact,  these  example  functions  may  provide  a  more  strenu¬ 
ous  test  of  the  forgoing  procedures  than  the  actual  functions  a{  T)  and  <(  T). 
Example  results  are  given. 

These  procedures  have  been  introduced  into  the  photon-electron  cas¬ 
cade  computer  program  EDPEC.  The  calculational  times  for  a  fixed 
number  of  cascades  have  been  reduced  to  about  50  percent  of  those  needed 
when  the  bisection  method  was  used  to  solve  the  integral  equations.  In  this 
comparison,  the  accuracy  criteria  for  each  case  was  set  to  be  much  higher 
than  the  overall  calculational  accuracy  of  EDPEC.  No  error  analysis  has 
been  conducted  to  determine  if  a  lesser  accuracy  (and  associated  decrease  in 
calculational  times)  can  be  used  without  degrading  the  overall  calculational 
accuracy  of  EDPEC. 

Smaller  improvements  in  the  calculational  efficiencies  may  still  be 
obtainable  by  devising  better  procedures  for  calculating  the  first  esti¬ 
mate  of  the  solution.  This  possibility  will  be  kept  in  mind  as  EDPEC  is 
used  in  future  studies. 
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